Mon. Not. R. Astron. Soc. 000, [Q-?? (2002) Printed 2 November 2010 (MN WT&t style file v2.2) 



Determining the extragalactic extinction law with SALT. 
II. Additional sample * 

i Ido Finkelman 1 ! , Noah Brosch 1 , Alexei Y. Kniazev 2,3 , Petri Vaisanen 2 ' 3 , 
C<l ; David A.H. Buckley 2,3 , Darragh O'Donoghue 2,3 , Amanda Gulbis 2,3 , Yas Hashimoto 4 , 
> ; Nicola Loaring 2,3 , Encarni Romero-Colmenero 2 ' 3 and Ramotholo Sefako 2 ' 3 

t O i 1 The Wise Observatory and the Raymond and Beverly Sackler School of Physics and Astronomy, the Faculty of Exact Sciences, 
|/ r < ■ Tel Aviv University, Tel Aviv 69978, Israel 

2 South African Astronomical Observatory, PO Box 9, 7935 Observatory, Cape Town, South Africa 
, 3 Southern African Large Telescope Foundation, PO Box 9, 7935 Observatory, Cape Town, South Africa 

4 Department of Earth Sciences, National Taiwan Normal University, No. 88, Sec. J^, Tingzhou Rd., Wenshan District, 
Taipei 11677, Taiwan R.O.C 



Accepted 2010 July 08. Received 2010 July 06; in original form 2010 June 14 



o 

u 

O ; ABSTRACT 

^ , We present new results from an on-going programme to study the dust extragalactic 

extinction law in E/S0 galaxies with dust lanes with the Southern African Large Tele- 
scope (SALT) during its performance-verification phase. The wavelength dependence 
of the dust extinction for seven galaxies is derived in six spectral bands ranging from 
C^l ■ the near-ultraviolet atmospheric cutoff to the near-infrared. The derivation of an ex- 

tinction law is performed by fitting model galaxies to the unextinguished parts of the 
image in each spectral band, and subtracting from these the actual images. We com- 
pare our results with the derived extinction law in the Galaxy and find them to run 
parallel to the Galactic extinction curve with a mean total-to-selective extinction value 
of Ry = 2.71 ± 0.43. We use total optical extinction values to estimate the dust mass 
for each galaxy, compare these with dust masses derived from IRAS measurements, 
. and find them to range from 10 4 to 1O 7 M0. We study the case of the well-known 

dust-lane galaxy NGC2685 for which HST/WFPC2 data is available to test the dust 
distribution on different scales. Our results imply a scale-free dust distribution across 
the dust lanes, at least within ~ 1" (~ 60 pc) regions. 

■ Key words: galaxies: elliptical and lenticular, cD; galaxies: ISM; dust, extinction. 
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INTRODUCTION 

Elliptical galaxies were considered for many years to be in- 
ert systems, devoid of interstellar matter. The study of ab- 
sorbing material in early-type galaxies gained interest with 
the identification of a class of galaxies characterized by an 
elliptical-like stellar body crossed by dark lanes or showing 
dark patches (Bertola h Galletta 1978; Bertola 1987). These 
are believed to be light extinguishing dust clouds. Much of 
our knowledge of the dust grain properties is derived from 

* based on observations made with the Southern African Large 
Telescope (SALT) 
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the analysis of the attenuation of starlight by the absorp- 
tion and scattering of light by grains in dusty structures. In 
particular, deriving the normalized amount of extinction as 
a function of wavelength, also known as the extinction law, 
is a useful tool to study the dust grain size and chemical 
composition. 

The extinction law in the Milky Way (MW) is usually 
studied by observing stars with similar spectral type and lu- 
minosity class but with different foreground extinction (Sav- 
age & Mathis 1979; Cardelli, Clayton & Mathis 1989). Since 
this method requires individual stars to be resolved, different 
methods must be invoked to study the extinction law even 
in nearby galaxies. In galaxies where the underlying light 
distribution across the stellar body varies smoothly, namely 
ellipticals and lenticulars, dust lanes and patches appear as 
local disturbances of the global brightness pattern. Mask- 
ing the dusty regions allows the extraction of a dust-free 
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model of the original underlying galaxy and the derivation 
of extinction at these regions. 

A significant number of E/SO galaxies have been inves- 
tigated with this method (see Goudfrooij et al. 1994; Patil 
et al. 2007; Finkelman et al. 2008; Finkelman et al. 2010). 
The overall result from these studies is that the extragalactic 
dust is amazingly similar to the dust in the MW, at least in 
the optical part of the spectrum. One possible interpretation 
of this result is that dust grains in different objects seem to 
have similar sizes and chemical composition. Such a similar- 
ity may have implications on the physical conditions and/or 
the origin of the dust at least for this class of galaxies. 

This paper is the second in our on-going program to 
observe dark- lane E/S0 galaxies with the Southern African 
Large Telescope (SALT; Buckley, Swart & Meiring 2006) 
and its CCD camera (SALTICAM, O'Donoghue et al. 2006). 
We present our analysis of seven dust- lane E/S0 galaxies, 
adding to the nine objects recently studied by Finkelman 
et al. (2008). Observing nearby dust-lane E/S0 with an (ef- 
fectively) nine-meter telescope has the specific advantage of 
achieving reasonable sensitivity even at the shortest spec- 
tral band used, despite the atmospheric extinction. The ex- 
tention of the optical study of the extinction law towards 
shorter wavelength allows the extraction of more informa- 
tion about the extragalactic extinction law close to the 
blue atmospheric cutoff where differences among extinction 
curves are better detected. 

This paper is organized as follows: §[T]gives a description 
of all the observations and data reduction; we present and 
analyze our results in §[2] discuss them in §[3j and summarize 
our conclusions in § [J 



1 OBSERVATIONS AND DATA REDUCTION 

CCD imaging observations of 12 dust-lane E/S0 galaxies 
were performed in service mode from February to May 2009 
with SALT during its Performance Verification (PV) phase. 
The observations were performed using the Johnson- Cousins 
B, V, R & I filters, as well as two short- wavelength interfer- 
ence filters, Ul and U2, with transmission bands peaking at 
340 nm (FWHM 35 nm) and 380 nm (FWHM 40 nm) respec- 
tively. The SALTICAM is a CCD mosaic of two 2048x4102 
pixel CCDs. The full SALTICAM image is about 9'.6 by 9'.6 
with pixels of ^0".28 (after on-chip binning by a factor of 
two), but the nominal science field is eight arcmin. 

The sample galaxies are listed in Table [1] with their co- 
ordinates, morphological classification and optical size taken 
from the RC3 catalog (de Vaucouleurs et al. 1992) or from 
HyperLEDA (Paturel et al. 2003). The objects were selected 
from catalogs of objects with dust lanes (Hawarden et al. 
1981; Ebneter & Balick 1985; Bertola 1987) and were chosen 
according to visibility conditions and availability of observ- 
ing nights, and not as part of a complete sample. 

The observations consisted of short-duration integra- 
tions while the filter wheel cycled between the six filters, re- 
sulting in six images per cycle. Typical exposure times were 
200 sec for the Ul & U2 bands, 30 sec for the B & V bands, 
and 15 sec for the R & I bands. They were determined to be 
sufficiently short to prevent image trailing when operating 
in the SALTICAM unguided mode (but with SALT track- 
ing) or overexposure of the brightest part of each galaxy. 



Due to the limited image quality of the telescope during 
the PV phase the workable field was limited to about four 
arcmin. The area nevertheless was sufficient to obtain min- 
imally overlapping multiple dithered images of the galaxies 
in the central regions of both CCDs. 

We created a sky flat-field from the set of exposures 
in each field for each filter by median-combining the scaled 
images without correcting for dithering. We then corrected 
the images for flat-field and created a count-frequency his- 
togram for each galaxy surroundings. We fitted a Gaussian 
distribution curve to the peak of the histogram, where the 
most frequent count value was assumed to represent the sky 
background. The images where the histogram significantly 
deviated from a Gaussian shape were rejected from the 
analysis. The reduced images were geometrically aligned by 
measuring centroids of several common stars in the galaxy 
frames and were then background-subtracted and combined 
to improve the S/N ratio. This alignment procedure involved 
IRAFlj tasks for translation and rotation of the images, so 
that a small amount of blurring was introduced, affecting 
the positional accuracy to a few hundredths of a pixel. 

While rsj 10—20 exposures were obtained for each galaxy 
in each filter, we selected only the images with the best im- 
age quality, and with seeing below 2.5". The final set of 
images for each filter therefore included only 3 — 5 images. 
Since the PSF shape of the stars in each image is often 
slightly distorted, we did not attempt to convolve the im- 
ages to the same FWHM before combining the images. This 
has no consequence for the results since the spatial scale of 
features we measure is determined by the image with the 
worst PSF size. The combined R-band contour map and B- 
R colour- index map for each galaxy are presented in Figures 
[T]and[2] respectively. 



2 ANALYSIS 

2.1 Dust extinction law 

The method for deriving the extragalaxtic extinction law 
in the dust lanes is described shortly here and we point 
the reader to Finkelman et al. (2008) for more details. The 
method assumes that elliptical galaxies have a smooth and 
symmetric light distribution with respect to their nuclei. It is 
therefore possible to fit these light distributions with a global 
brightness profile of a spheroid, e.g., a Sersic profile. We fit 
ellipses to the isophotes of the galaxy. Any disturbances of 
the global brightness profile of the underlying galaxy can 
be recognized by comparing the actual observed image with 
the model image. Suspicious dust features are expected to 
cause reddening due to the extinction of light by dust. The 
amount of extinction can be measured by computing how 
much light is missing in the suspected regions relative to 
the fitted smooth profile. 

The derivation of a dust-free model was performed in 
several steps. We fitted the galaxies with elliptical isophotes 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories (NOAO), which is operated by the Association of Uni- 
versities, Inc. (AURA) under co-operative agreement with the 
National Science Foundation 
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Figure 1. R-band contour maps. 
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Table 1. Global parameters for galaxies in our sample. 



Object 


RA 


DEC 


Morph. 


B° 


VHelio 


Size 




(J2000.0) 


(J2000.0) 


(RC3) 




(km/s) 


(arcmin) 


AM0052-321 


00:54:55 


-32:01:34 




14.90 


9607 


1.4x0.5 


NGC612 


01:33:19 


-36:26:38 


SOa 


14.04 


8995 


1.5x1.4 


AM0219-343 


02:21:19 


-34:19:08 


SO 


14.79 


6295 


1.4x0.9 


ESOH8-19 


04:19:00 


-58:15:27 


SO 


15.04 


1239 


0.9x0.8 


NGC1947 


05:26:48 


-63:45:36 


E-S0 


11.65 


1176 


3.5x3.2 


NGC3302 


10:35:47 


-32:21:31 


SO 


13.57 


3794 


1.9x1.3 


NGC3585 


11:13:17 


-26:45:18 


E 


10.82 


1434 


6.6x3.2 


NGC4753 


12:52:22 


-01:11:59 


SO-a 


10.85 


1242 


6.5x3.0 


NGC5266 


13:43:02 


-48:10:10 


SO 


12.06 


3707 


3.0x2.1 


NGC5363 


13:56:07 


+05:15:17 


SO-a 


11.10 


1136 


4.2x2.8 


AM1307-464 


13:10:33 


-46:59:27 


10 pec 


13.20 


400 


3.3x1.8 


AM1444-302 


14:47:27 


-30:38:43 


Sa 


14.71 


6649 


1.1x0.7 



Table 2. R\ extinction values. 



Object 


Rui 


RU2 


Rb 


Rv 


Rr 


Ri 


ESOH8-19 


4.84 ±0.60 


3.57 ±0.52 


3.50 ± 0.43 


2.50 ±0.31 


2.02 ±0.25 


1.46 ±0.18 


NGC1947 


4.60 ±0.31 


3.77 ±0.25 


3.53 ±0.23 


2.53 ±0.17 


1.90 ±0.17 


1.31 ±0.09 


NGC3302 


4.20 ± 0.40 


3.56 ±0.34 


3.52 ±0.33 


2.52 ±0.24 


1.03 ±0.11 


0.48 ±0.05 


NGC4753 


5.07 ±0.28 


4.33 ±0.24 


3.96 ±0.22 


2.96 ±0.17 


2.52 ±0.14 


1.90 ±0.11 


NGC5266 


5.35 ±0.82 


4.73 ±0.61 


4.41 ±0.57 


3.41 ± 0.44 


2.87 ±0.37 


2.36 ±0.30 


NGC5363 


4.82 ±0.36 


4.26 ±0.32 


3.98 ±0.29 


2.98 ±0.22 


2.62 ±0.19 


1.47 ±0.18 


AM0219-343 


5.32 ±0.55 


4.17 ±0.43 


3.09 ±0.28 


2.09 ±0.19 


1.75 ±0.16 


0.76 ±0.10 


MW galaxy 


4.90 


4.56 


4.10 


3.10 


2.32 


1.50 



using the ISOPHOTE package in IRAF while masking fore- 
ground stars and determining the center of the galaxy. The 
fit yielded for each isophote the ellipticity, position angle 
and mean isophotal intensity. We created a model unextin- 
guished galaxy from the parameters of the fitted elliptical 
isophotes, subtracted the observed image from the model 
image, and identified the dusty regions which we masked. 
We iterated the fitting procedure excluding all pixels with 
appreciable extinction (Av > 0.02, see e.g. Tran et al. 2001) 
and fixing the center if necessary. 

To measure the extinction we translated rectangular 
seeing-size boxes over the dust-occupied regions of each 
galaxy with no overlap. The amount of extinction in a par- 
ticular band was then calculated on a magnitude scale by 



-2.5 log 



I\,ob 



(i) 



while excluding the nuclear regions (radius ^ 5") to avoid 
seeing-related effects. No correction for Galactic extinction 
was made, assuming the Galactic extinction in the line of 
sight is uniform over each galaxy. 

By fitting a linear regression between the total extinc- 
tion A\ at wavelength A and the V-band extinction Av mea- 
sured along the dusty region in different bands we derived 
the extinction values for each galaxy. The best-fitting slopes 



were then used to calculate the 



Rx (e 



values for 



each galaxy and were compared with the Galactic values. 
The obtained R\ extinction values are listed in Table [2] 
along with the canonical extinction values taken from Sav- 



age & Mathis (1979) for comparison. The extinction curves 
are plotted in Figure [3] 

This method is limited to galaxies where the dust cov- 
ering factor is relatively small, and requires that the dusty 
regions are well-sampled. Therefore, for five galaxies where 
the apparent dust coverage is restricted to the inner parts of 
the galaxies or where the Sersic profile fits badly the galaxy 
we could not derive an extinction law as further described 
below. 

AM0052-321 shows a prominent dust lane along its ma- 
jor axis, but reveals a shell structure which prevents a proper 
modeling of the underlying galaxy. The optical image of the 
E/S0 galaxy AM1307-464 shows a patchy dust structure in 
its inner regions. However, the B-R colour image shows that 
this region is bluer than the underlying galaxy, implying 
that the colour variation is induced by a recent star for- 
mation episode and not by the presence of dust (see also 
Davidge 2007). Although AM1444-302 is classified Sa, it ap- 
pears in our images as an elliptical galaxy with a complex 
structure of dust lanes in its inner part. Since most of the 
inner region is obscured, and since the outer regions show 
shells, the underlying galaxy cannot be modelled. The B-R 
colour image of NGC3585 shows that the inner disk is red- 
der with respect to the bulge (see also Fisher et al. 1996), 
implying the presence of dust. Since the inner isophotes of 
NGC3585 are disky, suggesting the presence of a inner disk 
within the dominant bulge, we cannot model the underlying 
galaxy in that region by using our method. In addition, our 
images show that the dust lane across the bulge of NGC612 
actually follows a stellar disk at high inclination, resembling 
the well-known Sombrero galaxy. 
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Figure 3. Extinction curves for the program galaxies (solid lines) along with the canonical curve for the Galaxy (dashed lines) 
for comparison. Extinction values derived in Finkelman et al. (2010) are represented by open circles. The error bars are la 
errors. 




2.2 Dust mass estimation 



Assuming the extragalactic dust grains are uniform through- 
out each galaxy and are similar to those in the MW we 
use the Mathis, Rumpl & Nordsieck (1977) two-component 
model consisting of individual spherical silicate and graphite 
grains with an adequate mixture of sizes to calculate dust 
masses. The dust mass is then calculated by integrating the 
dust column density over the image areas S occupied by dust 
lanes, yielding: 

a + 

M d = SxH d = Sxl d x / -7ra 3 p d n(a)da. (2) 



where a is the grain radius; n (a) represents the size distri- 
bution of dust grains; a_ and a+ represent the lower and 
upper cutoffs of the size distribution, respectively; pa is the 
specific grain densities which is taken to be ~ 3 gr cm -3 for 
graphite and silicate grains (Draine & Lee 1984); Id repre- 
sents the dust column length along the line of sight; and tih 
is the hydrogen number density (see Finkelman et al. 2008). 

The value of the product Id x uh is inferred by measur- 
ing total extinction and calculating the extinction efficiency 
(Finkelman et al. 2008). The hidden assumption is that the 
dust forms a foreground screen for the galaxy and neglecting 
any intermix of dust and stars within the galaxy. The dust 
mass estimate from total extinction values provides there- 
fore only a lower limit to the real dust content of the host 
galaxies. 
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Alternatively, the dust mass can be estimated based on 
the observed far-IR dust emission assuming thermal equilib- 
rium. The relation between the dust mass and the observed 
far-IR flux density, I u , is 

Md = r pdD2 oJ£m (3) 

where pd and D are the specific grain mass density and dis- 
tance of the galaxy in Mpc, respectively; Q u and B u (Td) 
are the grain emissivity and the Planck function for the 
temperature Td at frequency zy, respectively (see Hilde- 
brand 1983). The dust grain temperature can be calculated 
from the IRAS flux densities at 60 /im and 100 pm using 

/ \0.4 

T d = ( J^H (Young et al. 1989). Hot circumstellar dust 
may also contribute to the 60 /im and 100 /im passbands, 
although this effect on the 60 /im and 100 pm flux densi- 
ties is only a few percent (Goudfrooij & de Jong 1995). The 
quantity apd/Qv can be calculated for silicate and graphite 
grains assuming it is independent of a for A ^> a (see Hilde- 
brand 1983; Finkelman et al. 2008). 

The calculated dust mass estimates represent lower lim- 
its, since IRAS was not sensitive to dust cooler than about 
20K, which emits mostly at A > 100 /im. Table [3] lists the 
estimated dust mass from the total optical extinction and 
the estimated dust mass based on IRAS flux densities taken 
from the catalog of Knapp et al. (1989) for bright early- type 
galaxies. We assume for the calculations a standard cosmol- 
ogy with Ho=73 km sec -1 Mpc -1 . 



2.3 Fine dust structure - the case of NGC2685 

The derivation of an extinction law is restricted by the im- 
age angular scale and therefore requires measuring a mean 
extinction value for each resolution element. However, the 
presence of dust-free "holes" or high-opacity "bricks" in the 
dust structure can significantly affect these measurements 
(see for example Keel & White 2001). 

To quantify the effects of dust on scales smaller than 
the - 0.1 - 1 kpc resolution provided by SALTICAM we 
use data from HST/WFPC2. Unfortunately HST did not 
observe any of our sample galaxies in the optical region; we 
therefore examine NGC2685, a well-known dust-lane galaxy 
for which HST data is available, as a representative example. 
NGC2685 was observed with WFPC2 in the B (F450W), V 
(F555W) and I (F814W) passbands. The total exposures 
were 2300 sec (B), 1000 sec (V) and 730 sec (I), split into 
halves. The individual exposures in each filter were com- 
bined using the CRREJ task in STSDAS for cosmic-ray 
recognition and rejection. 

While on ^kpc scales the dust in NGC2685 seems to be 
organized in a well-defined helical structure, HST/WFPC2 
images clearly reveal more complicated structures on much 
smaller scales. The ~ 0.1" PSF resolution provided by 
HST/WFPC2 allows deriving the extinction down to 2 x 2 
pixel boxes, and therefore studying the extinction on a spa- 
tial scale of ^ 5 pc. We plot in Figure 2] the relative dis- 
tribution of the B transmission values (=10 _As ) as mea- 
sured for 2 x 2 and 16 x 16 pixel boxes translated along the 
dust lane. Comparing the two relative distributions and the 
computed cumulative distributions hints at a rather smooth 
dust coverage and supports a scale-free structure on these 



scales. Furthermore, the B transmission values in the sub- 
structures vary from 0.2 to 0.9 and are consistent with low 
opacity within most of the obscured area (see for comparison 
Figure 6 in Keel & White 2001). 

The scale-free dust distribution can also be examined in 
terms of fractal analysis (see Mandelbrot 1982). A similar 
analysis of the interstellar dust structure in several nearby 
spiral galaxies revealed a hierarchical, fractal-like structure 
in a wide range of scales from ~ 0.1 pc to ~ 1 kpc (e.g., 
Westpfahl et al. 1999; Keel & White 2001; Sanchez & Al- 
faro 2008). Simply put, a fractal structure possesses details 
at every scale and a fractal dimension determines whether a 
system is homogeneous, and what fraction of space is filled 
(Combes 2000; Sanchez, Alfaro & Perez 2005). Astrophysical 
fractals are an approximation in a range between two limit- 
ing scales and are different from pure mathematical fractals 
that are infinite (Combes 2000). 

One method to test structural scaling or fractal be- 
havior is by constructing a scale-free relation between the 
perimeter of a contour and the scale length over which 
it is integrated or smoothed (Westpfahl et al. 1999; Keel 
& White 2001; Chappell & Scalo 2001). The box-counting 
procedure provides a simple perimeter-scale test (Westp- 
fahl et al. 1999). The scale- free relation, also known as the 
box-counting dimension, can be obtained by a power-law 
fit to the distribution of the number of boxes for different 
box sizes. Linear structures show a slope of —1, while frac- 
tal structures show a slope between —1 and —2. The box- 
counting dimension is the negative of this slope. 

Following the box-counting algorithm we used the ex- 
tinction maps to evaluate the length of single-total extinc- 
tion contours in the range Ab — 0.5 — 1.5 for smoothing 
lengths of 2, 4, 8 and 16 pixels, so that no resampling of 
pixel values affects the results. Figure \E\ plots on a log-log 
scale the perimeter calculated for each contour for each box 
size. The self-similar behavior of each contour is manifested 
as a straight line in this log- log plot, with slopes in the range 
~ ( — 1.1) — ( — 1-2) derived for the different contours. 

Determining the extinction law of a self-similar dust dis- 
tribution is not expected to be affected by the linear scale 
on which the extinction is measured. To verify this we de- 
rived the extinction law by measuring the extinguished light 
within 2x2 and 16 x 16 pixel boxes. The results are found to 
be consistent within the error bars and the mean extinction 
values are plotted in Figure [6] For comparison, we plot in 
Figure [6] also the extinction values measured by using the 
same method on SDSS images with ~ 1.2 arcsec resolution. 



3 DISCUSSION 
3.1 Dust extinction 

We derived extinction curves and extinction values for seven 
of our 12 sample galaxies, excluding AM0052-321, AM0131- 
365, NGC3585, AM1307-464 and AM1444-302. The mean 
extinction value for the galaxies, Rv — 2.71 ±0.43, matches 
the rsj 2.80 value derived in previous studies (e.g., Patil et 
al. 2007; Finkelman et al. 2008) for E/S0 galaxies with well- 
settled dust lanes and is slightly smaller than that in the 
MW, but similar within the error bars. These Rv values are 
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Table 3. Dust mass. 



Object IRAS flux (mjy) Log^)^ Log 





60/im 


100/xm 






ESOH8-19 


733 ± 51 


1391 ± 84 


4.70 ±0.08 


4.44 ±0.02 


NGC1947 


1052 ± 53 


4136 ± 248 


5.69 ±0.09 


5.01 ±0.01 


NGC3302 


1876 ± 206 


5401 ± 648 


6.56 ±0.17 


5.55 ±0.02 


NGC4753 


2438 ± 293 


9008 ± 1171 


6.02 ±0.20 


5.20 ±0.01 


NGC5266 


1220 db 98 


4365 ± 349 


6.63 ±0.12 


5.04 ±0.03 


NGC5363 


1693 ± 169 


5150 ± 567 


5.54 ±0.16 


4.86 ±0.02 


AM0219-343* 


0±42 


0±210 


< 8.2 


6.00 ±0.02 



*Upper limits assuming T = 20K 



Figure 4. Distribution of B-band transmission by area. The up- 
per panel shows the differential distribution of dusty regions with 
different optical depth values while translating a 2 x 2 pixel box 
(~5 pc; solid line) and a 16 x 16 pixel box (~40 pc; dashed line) 
along the dust lanes. The data suggest that most regions have B 
transmission from 0.2 to 0.9. The lower panel refers to the cu- 
mulative distribution of the data presented in the upper panel. 



Figure 5. Box counting analysis of the relation between length 
(region perimeter) and smoothing length for a 256 x 128 pixel 2 re- 
gion in NGC2685. Single contours were evaluated for total extinc- 
tion values from 0.5 to 1.5 represented in the plot by a different 
symbol with the lines marked by the extinction value. Scale-free 
(fractal) behavior is manifested as straight lines in this log-log 
plot. 




0.4 0.6 0.8 

B transmission 



1.2 




0.4 0.6 0.8 
Log(box size) 

Galactic values are generally linked with a variation of the 
characteristic grain size responsible for the optical extinc- 
tion in such galaxies with respect to dust in the MW (see 
Finkelman et al. 2008 and references therein). However, it 
should first be established that the measured extinction law 
represents the "true" extinction law, as discussed below. 



also within the range of values derived by Goudfrooij et al. 
(1994), 2.1 to 3.3, which implies a characteristic dust grain 
size up to 50% smaller than that of standard Galactic dust 
grains (see Finkelman et al. 2008). The lowest extinction 
value measured in a galaxy in our sample, Rv = 2.09 ± 
0.19, is similar to the lowest extinction values obtained by 
Goudfrooij et al. (1994), Patil et al. (2007) and Finkelman 
et al. (2008). These results further support Goudfrooij et al. 
(1994) conclusion that E/S0 galaxies with smooth and well- 
defined dust lanes exhibit smaller Rv values (see also Patil 
et al. 2007; Finkelman et al. 2010). 

We show that although the derived extragalactic extinc- 
tion curves are close to the standard Galactic value, most 
of them lie above or below the Galactic extinction curve 
and parallel to it. Deviations of extinction values from the 



3.2 Effect of young stars 

Goudfrooij et al. (1994) discussed the effect star formation 
and light emission from young massive stars on the derived 
extinction law. Finkelman et al. (2010) studied a sample of 
galaxies where the dusty structure is generally morphologi- 
cally associated with ionized gas; this hints at the presence of 
a small fraction of young stars and of a recent star formation 
episode. Since the light of the underlying galaxy is extrapo- 
lated from the unobscured regions, and since the presence of 
a young stellar population in the dusty locations is expected 
to be bluer than the old stellar population, the extinction 
might be underestimated, especially at the shorter wave- 
lengths. A uniform young stellar population in the dust lane 
will have the same effect over all obscured regions and there- 
fore should not affect the slope of the A\ vs. Ay relation, 
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Figure 6. Extinction curve for NGC2685 from HST (upper 
curve) and SDSS (lower curve) data. The extinction values (plus 
signs) are presented along with the canonical curve for the Galaxy 
(dashed lines) for comparison. The error bars are la errors. 




r 1 (Jim- 1 ) 



but only its intercepts with the axes. However, with a more 
complex stellar structure the real extinction curve may be 
more concave than the measured one and show larger ex- 
tinction values. 

The effect of forward scattering in dusty regions on the 
observed extinction curve can be neglected as well. This ef- 
fect is mainly dominant at the edges of disks and leads to 
apparently lower extinction values in the blue and to a rel- 
atively flat extinction curve with respect to the Galactic 
curve, whereas such flattening is not evident in any of our 
sample galaxies (see Finkelman et al. 2008). 



3.3 Effect of resolution and scales of dust 
distribution 

The PSF of the SALT observations is typically 2.5". The 
extinction is therefore measured with an assumption that 
the dust structure is homogeneous and smooth within each 
resolution element. Furthermore, the hidden assumption is 
that A\ = 1 .086 t \ is valid for a given wavelength only 
as long as the dust is a foreground screen, while the op- 
tical depth r\ may vary with column density and column 
length. Finkelman et al. (2008) studied the influence of the 
dust location on the observed extinction law for the case 
of varying optical depths concluding that if the dusty lo- 
cations are optically-thin, i.e., tv — 1 as typical for our 
sample galaxies, the deviation from the "real" linear rela- 
tion may not be significant. However, tracing the possible 
variations of the dust column density and geometry within 



the dust lane in a certain galaxy requires better spatial res- 
olution than provided by SALTICAM. With no available 
high-resolution optical data of our sample galaxies, we used 
archival HST/WFPC2 data of NGC2685, a representative 
E/S0 galaxy with well-defined dust lanes, to test the Finkel- 
man et al. (2008) hypothesis. In particular, we investigated 
the presence of highly obscured regions through the search 
of greyish extinction. Our analysis shows that the dust lanes 
do not show transparent "holes" or opaque "bricks" , but are 
rather smooth structures down to scales as small as 5 — 10 pc. 
Figure [3 further demonstrates the foreground screen model 
validity for the case of our galaxies. 

The dust structure in NGC2685 is characterized by a 
fragmented and filamentary geometry. The complexity of 
structures of interstellar clouds are often studied and charac- 
terized by fractal analysis (see Khalil et al. 2006; Sanchez & 
Alfaro 2008 and references therein). The limited resolution, 
finite-size effects and use of different analysis algorithms are 
known to introduce biases on the estimate of the fractal di- 
mension. Using the box-counting technique, previous studies 
of molecular and dust clouds found that the fractal dimen- 
sion value is almost always between 1.2 and 1.5. 

Although a significant part of fractal analysis studies is 
dedicated to the distribution of gas and dust in the Galaxy 
(e.g., Khalil et al. 2006), a growing number of studies also fo- 
cus on the HI (e.g., Westpfaul et al. 1999), HII (e.g., Sanchez 
& Alfaro 2008) and the dust component (e.g., Keel & White 
2001; Elmegreen, Elmegreen & Leitner 2003) in nearby spi- 
ral galaxies. Therefore, the general result of a more or less 
constant fractal dimension might reflect a universal value 
that is related with the underlying physical processes that 
drive the growth and subsequent evolution of the interstel- 
lar clouds (e.g., Elmegreen & Falgarone 1996; Elmegreen 
1997; Combes 2000; Chappell & Scalo 2001). The similar- 
ity in fractal structure of both dust and molecular clouds 
also support an hypothesis that the dust is related with the 
process of molecule formation in cold clouds (Datta 2004). 

Examining the extinction map of NGC2685 we find the 
contours at constant extinction to be individually scale-free 
for sampling lengths ^ 2 pixels (^ 5 pc). We measured for 
NGC2685 a box-counting dimension between 1.1 and 1.2, 
which is only marginally consistent with the fractal dimen- 
sion measured for the ISM in previous studies, and is closer 
to the fractal dimension of linear and euclidean structures. 
Therefore, it is not clear whether the fractal behavior of 
ISM in dust-lane E/S0 galaxies follows the universal struc- 
tural evolution of ISM observed in the MW and nearby spiral 
galaxies or otherwise indicates a different evolution. Deter- 
mining this with confidence requires measuring with HST 
the fractal dimension of the dust structures in a larger sam- 
ple of dust-lane E/SOs. 

3.4 Dust mass 

The dust mass for each sample galaxy was evaluated inde- 
pendently in two ways, based on the optical extinction or 
from far-IR observations. The estimated masses lie in the 
range ~ 10 4 to 1O 7 M0, in good agreement with previous es- 
timates for early-type galaxies with dark lanes (see Finkel- 
man et al. 2008 and references therein). Since we could not 
measure the entire dusty region with the optical method, 
and since we assumed that the dust does not intermix with 
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Figure 7. Left panel: Extinction values measured in NGC1947 from SALT data. Right panel: Extinction values calculated from 
HST/WFPC2 by measuring the intensity in 2x2 (grey plus signs) and 16x16 pixel boxes (black plus signs) translated across the dusty 
regions of NGC2685. The solid lines represent the model best-fit assuming dust is a foreground screen. 
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the stars within the galaxy, this method provides only a 
lower limit to the true dust content of the host galaxies that 
should be confirmed by other methods. Using IRAS flux 
densities to trace the dust far-IR emission reveals that the 
dust masses are up to an order of magnitude higher than 
measured from the optical extinction, implying the presence 
of a more massive diffuse dust component (Goudfrooij & de 
Jong 1995; Patil et al. 2007; Finkelman et al. 2008). 



4 CONCLUSIONS 

This is the second paper presenting new results from an 
on-going program to study the extragalactic extinction law 
with SALT down to the ultraviolet atmospheric cutoff (see 
Finkelman et al. 2008). 

We derived the ratio between the total V band extinc- 
tion and the selective B and V extinction Rv for seven 
out of our 12 sample galaxies, and obtained an avergage 
of 2.71 =b 0.43, slightly smaller than that found by Finkel- 
man et al. (2008) and than the Galactic value but similar 
within the error bars. Extinction curves derived for the sam- 
ple galaxies run parallel to the canonical Galactic curve but 
do not necessarily overlap. We found that deviations from 
the standard Galactic value are unlikely to be affected by 
the image resolution and scales of the dust distribution, and 
that the effect of intermix of stars and dust, which needs 
to be considered when measuring the extinction in highly- 
inclined disk galaxies, is negligible for the case of dust lanes 
in E/S0 galaxies. We conclude that dust in the extragalac- 
tic environment has similar properties to the Galactic grains, 
and that any deviation of extinction value from the Galactic 
one can be accounted for by a variation of the characteristic 
grain size responsible for the optical extinction. 

The dust content in the dust lanes derived by the optical 
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extinction method for the sample galaxies is in the range 
~ 10 4 to 1O 7 M0, whereas the diffuse dust content might be 
larger by an order of magnitude. These results are fully in 
agreement with those of Finkelman et al. (2008). 
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